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We present a model of pressure effects of a two-band superconductor based on a Ginzburg-Landau 
free energy witfi two order parameters. The parameters of the theory are pressure as well as tem- 
perature dependent. New pressure effects emerge as a result of the competition between the two 
bands. The theory then is applied to MgB2. We identify two possible scenaria regarding the fate 
of the two G subbands under pressure, depending on whether or not both subbands are above the 
Fermi energy at ambient pressure. The splitting of the two subbands is probably caused by the E2g 
distortion. If only one subband is above the Fermi energy at ambient pressure (scenario I), applica- 
tion of pressure diminishes the splitting and it is possible that the lower subband participates in the 
superconductivity. The corresponding crossover pressure and Gritneisen parameter are estimated. 
In the second scenario both bands start above the Fermi energy and they move below it, either 
by pressure or via the substitution of Mg by Al. In both scenaria, the possibility of electronical 
topological transition is emphasized. Experimental signatures of both scenaria are presented and 
existing experiments are discussed in the light of the different physical pictures. 
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I. INTRODUCTION 



The discovery of superconductivity in the material 
MgB2 [1] initiated intensive recent theoretical and ex- 
perimental interest. The possibility of a high critical 
temperature in a class of materials which are chemi- 
cally much simpler than the high-Tc cuprates and the 
occurence of large critical current densities pose some in- 
teresting new questions in the research of superconduc- 
tivity. MgB2 is a type II superconductor in the clean limit 
[2] . The crystal belongs to the space group P6 jmmni or 
AZi?2-structure where borons are packed in honeycomb 
layers alternating with hexagonal layers of magnesium 
ions. The ions Mg^+ are positioned above the centers of 
hexagons formed by boron sites and donate their elec- 
trons to the boron planes. The electronic structure is or- 
ganized by the narrow energy bands with near two-fold 
degenerate (cr-electrons) and the wide-band (7r-electrons) 
[3-5] . Without any of the lattice strain, the u dispersion 
relations are slightly splitted due to the two boron atoms 
per unit cell (the electronic analog of Davydov splitting 
for phonons). The a portions of the Fermi surface (FS) 
consist of coaxial cylinders along the F - A symmetry di- 
rection of the Brillouin zone (BZ), whereas the tt- bands 
are strongly dispersive. 

The excited phonon modes in MgB2 present a sharp 
cut-off at about 100 meV. The optical modes (Big, E^g., 
and are practically non-dispersive, along the 

T — A direction. Various Raman data [6-8] show a small 
spread of the i?2g frequency around 74.5 meV in differ- 
ent ceramic samples of MgB2. The strong deformation 
potential A of the in- plane E^g mode [9,10] causes a sig- 
nificant energy splitting of the o band around 1.5 eV, 



lifting its two-fold degeneracy at the F-point of the BZ. 
In general, measurements of the effect of pressure on the 
electronic structure are based on the Raman technique, 
therefore only the influence on the Eig mode can be 
traced. High pressure experiments up to 15 GPa [6] have 
revealed a large increase of the out of phase E^g phonon 
mode. This leads to a suppression of the displacement 
of the boron atoms as well as of the deformation poten- 
tial A = B^u^ + B^u^ where u is the displacement [10] 
which will be used below. Moreover, experiments under 
pressure up to 40 GPa do not show any structural phase 
transition [6,11-13]. Due to the anisotropy the hydro- 
static pressure decreases the inter-plane distance more 
than the in-plane distance between adjacent borons. 

The interpretation of experimental data from different 
spectroscopic methods [14-16], suggests the presence of 
two different superconducting gaps. The specific heat 
behavior [17], the low isotope effect [18], pressure effects 
[19] and penetration depth data [20] provide evidence 
for a complicated superconducting order parameter in 
MgB2. Recent NMR measurements [21] of as well 
as Hall measurements [22] are consistent with a a-band 
driven superconductivity in MgB2 where the tt band par- 
ticipates due to interband scattering. The use of two or- 
der parameters is justified by recent experimental results 
on scanning tunneling microscopy [23]. 

Some microscopic theories [9,24] are based on a band 
or cr — TT band scenaria for superconductivity whereas 
others [25] are concentrated on a tt band superconduc- 
tivity. A consensus is formed though, that the driving 
band for the superconductivity is the cr one. The cal- 
culated spectral functions [26] and the analysis of the 
reflectance measurements [27] show the possibility of dif- 
ferent superconducting mechanisms beyond the conven- 
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tional electron-phonon Bardccn-Coopcr-Schrieffer (BCS) 
pairing. The complicated nature of superconductivity 
in MgB2 does not allow at present to accept the con- 
ventional superconducting mechanism. Therefore, it is 
appropriate to construct a phenomenological theory. 

Since MgB2 presents an example of a two-band super- 
conductor [28] may serve as a paradigm to investigate 
new effects. The purpose of this article is to present 
a phenomenological Ginzburg-Landau (GL) theory for a 
two-band superconductor appropriate to MgB2 based on 
two order parameters and, in that framework, to provide 
a physical picture of pressure effects on MgB2. The study 
of pressure effects is one way to investigate (i) the differ- 
ent topology of the two bands and (ii) the question of 
the participation of the two a subbands in the supercon- 
ductivity. A brief account of a part of this work was pre- 
sented elsewhere [29] . In the present article missing terms 
in the GL functional are restored, non-adiabatic effects 



are explained using a different approach and the effects 
of pressure are discussed in more detail and compared 
with very recent experiments. The work is organized as 
follows. In section II the GL description is presented, in 
Section III the pressure effects, in connection with the a 
band, are discussed and a general discussion with some 
concluding remarks are given in Section IV. 



II. GINZBURG-LANDAU DESCRIPTION 

We introduce the GL free energy functional with two 
order parameters, appropriate for MgB2, assuming that 
both the order parameters belong to the Fi representa- 
tion of the point group of MgB2 crystal. The two order 
parameters are labeled by the corresponding band (cr or 
tt) of MgB2, without loss of generality : 



+ 7i(n,v.n:^; + n^V'an;^: + c.c.) + /3i^<.|2|V.l' +^{yxAf}, (i) 



where fl = —ihV — 2e/cA is the momentum operator, A 
is the vector potential and ttg-.Tr = (^^a ti{T — T^^. ^). It is 
possible to take into account the largo anisotropy of the 
two order parameters (the a is almost two-dimensional) 
by rescaling of the axis according to the effective masses 
(directional dependence effective mass). In the present 
work, since we are focusing exclusively on pressure ef- 
fects, the derivative terms are not important. 

In Eq. (1) we have used the fact that the mixing r- 
term favors the coupling of linear combination of the two 
order parameters with phase difference if r < or tt if 
r > between them [30], therefore a term + c.c. 

is incorporated into the /3-term of the free energy [31]. 
If r = then Eq. (1) is the free energy for two bands 
without Josophson coupling between them. The quartic 
/j-term is the only one which mixes the two gaps and 
is unimportant in this case. The onset of the supercon- 
ducting state in one band does not imply the onset in 
the other. This corresponds to the case of Vsd = in the 
two-band BCS treatment [32]. If r ^ 0, which is the case 
of MgB2, then the pair transfering term is present in the 
GL functional and it means that the onset of supercon- 
ductivity in one band implies automatically the appear- 
ance of superconductivity in the other. There is a single 
observed Tc which is a function of the bare ones T°^^. 
In MgB2 the compression due to pressure is anisotropic 
[6,33,34]. According to Ref. [6] the compressibility along 
the c-axis is almost twice larger than that of the plane 
compressibility. Therefore application of uniaxial pres- 
sure on single crystals will affect differently the two gaps. 
Following Ozaki's formulation [35,36], we may add to the 



GL functional, Eq. (1), the term which couples the order 

parameters, in second order, with the strain tensor e to 
first order, having already specified the symmetry of the 
order parameters : 

Fstrain = " Ci(Ti)[5(exx + eyy) + ezz]\ipcr\'^ 

- Cz{Vr)[5'{exx + Cyy) + e,,]{ilj*7p^ + c.c), (2) 

where Ci_2.3(ri) arc coupling constants and S, S' are 
given in terms of the elastic constants of the material. 
This will load to a change of the bare critical tempera- 
tures T^^ ^ and consequently of the actual critical tem- 
perature Tc- The physical reasoning behind this is that 
the material will change in such a way as to gain con- 
densation energy by enhancing the density of states in 
the direction where the superconducting gap is larger. 
The term proportional to V'^V'ti- + cc. in the Fstrain will 
renormalize the mixing coupling r to f. We model the 
corresponding differences by taking a linear dependence 
of the two bare Tc's with pressure Tco-.t = T^a n ^ Va-,TrP , 
where ricr,Tr = \dT^^j^{p = 0)/dp\ and we will discuss the 
validity of the assumption for the MgB2 later. 

Analyzing the equations which result from the mini- 
mization of the GL free energy we get : 

a^aa=r'^. (3) 

This gives the pressure dependence of the superconduct- 
ing critical temperature: 



2 



+ l{[TL-T^. + iv.-V.)p? + a'V/^ (4) 

where = 4f^/(a°a°). Deviations from a straight 
hne at moderate values of pressure can be attributed to 
the two bands. From the above formula the inequality 
dTc/dp < is always true as long as the renormalized 
mixing parameter f is real, as a consequence of the ini- 
tial assumption on the pressure dependence of the bare 
critical temperatures. In the case of MgB2 we can safely 
consider that ry^ > rja- 

III. NONADIABATIC EFFECTS 

We now wish to address an issue specific to MgB2 
which can be also realized in other multiband supercon- 
ductors. The physical situation we would like to question 
is the splitting of the two a subbands at ambient pressure 
and their participation in superconductivity with the in- 
crease of pressure. Due to contradictory experimental as 
well as theoretical results we proceed by distinguishing 
two different scenaria. The first one addresses the split- 
ting of the two sigma subbands at ambient pressure as 
shown by first princible calculations [9,10,37] when one of 
the subbands is below the Fermi energy Ep. The second 
scenario addresses the opposite situation where both the 
a subbands are above Ep s± ambient pressure. 

A. Scenario I : one a sub-band above Ef at ambient 
pressure. 

In particular in [37] the splitting of the two subbands 
was studied in detail and compared with similar situation 
in AIB2. The conclusion is that due to the E2g phonon 
mode, the two cr bands split nearby T point and the lower 
band completely sinks below the Fermi energy. This is 
the case for MgB2 and also for the heavily hole-doped 
graphite. Moreover, experimental results [33,38] showed 
a kink in the superconducting critical temperature as a 
function of pressure at approximately 6-8 GPa and also 
a kink in the volume dependence of Tc for Mg^°B2 at 
around 20 GPa and Mg"B2 at around 15 GPa [7]. 

The physics we discuss here, is the change in the elec- 
tronic properties of the material under pressure and its 
influence on the superconducting state. The band which 
is below the Fermi level at ambient pressure [9] is possi- 
ble to overcome the energy difference and to get above 
the Fermi level at a certain value of pressure (termed as 
crossover pressure), restoring the degeneracy of the two 
a bands at point T. This is a non-adiabatic effect as dis- 
cussed in [37] . To understand and illustrate this effect on 
the superconducting state we need to consider the pres- 
sure dependence of the coefficient as usual [39] with 



the modification due to the particular physical situation 
by writing: 

a„ = a° (T - + r,„p) + a^p - p,)Q{p - p,), 

(5) 

where Q{x) is the Heaviside step function, and a J. are 
positive constants. This choice reflects the fact that at 
a certain crossover pressure Pc the second band starts to 
participate in superconductivity as well and that initially 
it is an increasing function of pressure. We obtained the 
above formula by requiring the continuity of the coeffi- 
cient a„ at Pc. 

After the energy difference between the two a bands 
almost disappears, then they both follow the same re- 
duction of the bare Tc's with pressure. The approximate 
value of the crossover pressiire can be estimated as fol- 
lows. The energy differc;nc;c! between the two subbands is 
approximated by 5E ~ (1 — n) a/A, where the fraction 
of superconducting electrons is 1 — n ~ 0.03 [3], i.e. the 
carrier density per boron atom in MgB2. The approxi- 
mate value of Pc is half the value of the pressure which 
suppresses the deformation potential A. The deforma- 
tion potential in turn, can be estimated by the expression 
Pc^ ^ (1 — n) \/A where Vl ^ 30A^ is a unit cell volume 
of Mgi?2 and A = 0.04ey^ is the deformation potential 
for a boron displacement u ^ O.OSA [10]. Using these 
parameters, we get a crossover pressure of Pc ~ 15GPa. 
This estimate shows, that a realistic applied pressure in- 
fluences drastically the electronic structure and the FS 
topology, restoring the degeneracy of the two subbands 
at the F-point which are initially splitted. Superconduc- 
tivity in the second subband may occur at lower values 
of the estimated Pc due to the fact that both subbands 
are affected. The elimination of the energy difference 
around the Fermi energy will also occur at lower values 
due to the corrugation of the a portions of the FS, the al- 
ready existing strains in the material and the anisotropic 
compressibility. These considerations make the above es- 
timate an upper limit of the crossover pressure pc- The 
degree to which shear stresses of the sample and the sur- 
rounding fluid imder pressure affect the data of pressure 
measurements is still under investigation [19]. 

The microscopic nature of this crossover is caused by 
the change of the FS. If the topology of the FS changes 
under pressure (this subsection) or due to substitutions 
in the composition (next subsection), a van Hove sin- 
gularity at critical energy Ec in the electronic density 
of states p{E) is manifested in an electronic topological 
transition of the 5/2 kind [40]. Since dTc/dp is essentially 
proportional to the derivative dp{Ep)/dEp, it develops 
singularities oc l/\/Ep — Ec, provided that the strength 
of carrier attraction varies slightly with pressure [41]. 

In Fig. 1 using a model calculation, we illustrate the 
expected behavior of Tc and the form of the order param- 
eter as a function of pressure at fixed temperature and. 
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schematically, the kink of at pc due to the change of the 
slope from Eq. (5). The predicted kink close to Pc can be 
detected directly in a penetration depth experiment un- 
der pressure. The chosen parameters are such that they 
respect the relation tIt^ > Vo-j the fact that the supercon- 
ducting density coming from the a band is higher than 
the one coming from the tt band (a^ > ctTv) and that the 
driving band for superconductivity is a (T^ > T^). 
The contribution on the Griineisen parameter : 

G^S^ = ^^, (6) 
dp Lu dp 

where B is the bulk modulus and uj is the phonon fre- 
quency, can be approximately obtained if we observe 
that SE ~ hujE2g then G ~ B/2pc. The bulk modu- 
lus is measured to be B~ 114 GPa [42] and our estimate 
for the Griineisen parameter is G~ 3.8. Experimentally 
G= 2.9 ± 0.3 as reported in [6] for the measured Raman 
active E2g phonon mode. However, as indicated in Ref. 
[6] , for anisotropic crystals it would be more appropriate 
to scale the frequency shift in the Griineisen parame- 
ter with the variation of the lattice constant a such as 
G = dlnuj/Sdlna [43]. Then, the Griineisen parameter 
for the MgB2 takes the value 3.9 ± 0.4. There is an ex- 
cellent agreement with the above estimate and it justifies 
our approach. 




P^Pc 

FIG. 1. The form of the order parameter as a function of 
pressure, where the kink at Pc is shown. Inset: the behavior 
of the critical temperature Tc(p). The chosen parameters 
are (in reduced units): T=0.8 T°^, T°„=0.7 T°^, 7]^= 0.2 
TiJPc, 0.5 T°Jp„ a«=2, a«=l, al= 0.1. 

B. Scenario II : both a sub-bands above Ef at 
ambient pressure. 

Some recent experimental and theoretical studies [44] 
lend support to a second scenario. More specifically, the 



substitution of Mg by the higher valent Al was used to 
provide the necessary change of the Fermi energy via elec- 
tron concentration and the lattice parameters in the com- 
pound Mgi_j;Ala;B2. This is in essence equivalent to a 
change of Fermi surface topology as in the pressure ex- 
periments. As it was shown in [44] by increasing the Al 
content (equivalently by increasing the pressure), a topo- 
logical crossover occurs at some point and the negative 
slope of Tc as a function of Al content x decreases further 
giving a different picture from that of Fig. 1. 

To model this in the language of the GL model, one ba- 
sically has to substitute the pressure p and the crossover 
pressure pc with the Al concentration x and and the 
crossover concentration Xc respectively, requiring again 
continuity at Xc- 

(^cr = a° ,.{T - T^^ + i^^^^x) + al^^{x - Xc)Q{x - Xc) 

(7) 

where the extra substrcipt x denotes that the numerical 
coefficients are different from Eq. (5). The critical tem- 
perature as a function of the Al content is illustrated in 
Fig. 2. The physics is similar to the topological transition 
described in scenario I where now instead of pressure, the 
changing parameter is the composition. It is a more con- 
trolled way to change the topology of the FS. In fact, an 
abrupt topological change in the cr-band Fermi surface 
was found at a; = 0.3 in Ref. [44]. When the a bands are 
filled in Mgi_a:Al3;B2 at x « 0.6, superconductivity diss- 
apears. It is worth noticing that in this compound, the 
impurity scattering broadens the van Hove singularities 
at the saddle points Ec of the FS and smears the change 
of slope of Tc at Ec. Therefore the kink of Tc in Fig. 2 
will be less pronounced (see also the experiment in Ref. 
[44]). 

We emphasize here the need of de Haas-van Alphen 
(dHvA) data for Mgi_j;Ala;B2, especially close to x=0.3. 
The pressure derivative of the superconducting critical 
temperature as a function of the concentration dTc{x)/dp 
can be also conclusive since it can clarify the correlation 
between pressure, Al concentration and the kink at Xc- 




FIG. 2. Schematically the form of the superconducting 
critical temperature Tc(x) in the second scenario, x is the 
concentration of Al in a Mgi_a;Ala;B2-like system. 
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IV. DISCUSSION 

Interestingly enough, there is an obvious discrepancy 
from the two pictures of the previous section, which calls 
for more experimental data. Observation of quantum 
oscillations in dHvA experiments provides the informa- 
tion about the bands, the effective mass of the carriers 
and the shape of the Fermi surface in MgB2. In Ref. 
[45], only the smaller cylindrical Fermi surface along the 
c-axis was observed (with a hole band mass -0.251 me, 
where me is the electron mass). The same authors re- 
ported the observation of only the larger of the two a 
tube-like Fermi surfaces (with a band mass -0.553 me) 
in Ref. [46]. Although the estimated parameters of the 
MgB2 are in agreement with the LDA band calculations 
[47] future dHvA experiments should clarify the detailed 
observation of the a portions of the FS. Also note that 
angle resolved photoemission spectroscopy experiments 
[48] detected one band along the T — K symmetry line of 
the BZ. 

There are a few points in order to discuss further with 
respect to applicability of the theory developed here. 
Early pressure experiments [13] demonstrated the over- 
all decrease of Tc with a pressure increase which was at- 
tributed to the loss of holes. In two of the samples of 
[13] there is a linear dependence of Tc on pressure and in 
two others a weak quadratic dependence. We stress that 
the samples were polycrystalline and the experiment is 
effectively under hydrostatic pressure. Also the degree 
of nonstoichiometry was not known. The almost linear 
dependence for a wide range of pressures, makes the GL 
functional as presented, valid for MgB2. 

Experiments on single crystals will be able to clarify 
the effects which arc described. We do not attempt at 
the moment any actual detailed fitting of experimental 
data, except from the illustrative fitting in section III, 
because there is no experimental consensus on the differ- 
ent values of key parameters of the theory {e.g., there is a 
wide range of published data on the value oidTc/dp [19]). 
More detailed Raman data with optical reflectivity, spe- 
cific heat and thermal conductivity measurements under 
pressure in the superconducting phase are needed as well 
as a penetration depth experiment which can reveal the 
pressure dependence of the superfluid density. Also the 
experiments suggested in scenaria I and II will be cru- 
cial in the question regarding the topological transition, 
where Ep can be varied experimentally either by exter- 
nal pressure (Sec. IIIA) or by alloying (Sec.IIIB). Recent 
Raman spectra at high pressure [38] reveal a reduction of 
the slope of the pressure-induced frequency shift by about 
a factor of two, at about 18 GPa which supports the sug- 
gestion that MgB2 may undergo a pressure-induced topo- 
logical electronic transition. In connection with the ap- 
pearance of superconductivity under pressure, we stress 
the high pressure superconducting phase of CaSi2 which 



superconducts above an applied pressure of 12 GPa with 
Tc ~ UK [49]. 

One proposed experiment which can potentially show 
the different role of the two bands (cr and tt) is the de- 
tection of the Leggett mode or the internal Josephson 
current [50]. Moreover an upward curvature in Hc2 is 
also known to be a signature of multicomponent super- 
conductivity. 

In summary, we present an analysis of pressure effects, 
within the GL theory, of a two-band superconductor and 
apply it to MgB2 . We make predictions for non-adiabatic 
effects and discuss different experiments from which cru- 
cial information can be extracted on the physics of the 
two participating bands cr and tt as well as the more del- 
icate questions on the two a subbands. 
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